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Abstract 

Measurement of mixing-induced CP violation in i?^ decays is of prime importance 
in probing new physics. So far only the channel i?^ — )• J/'ip(j) has been used. Here we 
report on a measurement using an LHCb data sample of 0.41 fb~^, in the CP odd 
eigenstate J/^/^/o(980), where /o(980) — ^ vr+vr". A time dependent fit of the data 
with the lifetime and the difference in widths of the heavy and light eigenstates 
constrained to the values obtained from i?^ — )• J/tp(j) yields a value of the CP 
violating phase of —0.44 it 0.44 it 0.02 rad, consistent with the Standard Model 
expectation. 
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1 Introduction 

An important goal of heavy flavour experiments is to measure the mixing-induced CP 
violation phase in 5° decays, 0^- As this phase is predicted to be small in the Standard 
Model (SM) P, new physics can induce large changes [2]. Here we use the decay mode 
Bl — >■ J /ijjfQi^SQ). If only the dominant decay diagrams shown in Fig. [l] contribute, then 
the value of 0^ using — J /ipf o{980) is the same as that measured using i?° — )■ J/ipcp 
decay. 




or K K" 



Figure 1: Dominant decay diagrams for 5^ — )■ J/?/'/o(980) or J/ipcf) decays. 

Motivated by a prediction in Ref. [3], LHCb searched for and made the first observation 
of 5° — )■ J/V'/olQSO) decays [3] that was subsequently confirmed by other experiments 
[5l |6]. Time dependent CP violation can be measured without an angular analysis, as the 
final state is a CP eigenstate. From now on /q will stand only for /o(980). 

In the Standard Model, in terms of CKM matrix elements, 0^ = —2 arg ^^yf • The 

equations below are written assuming that there is only one decay amplitude, ignoring 
possible small contributions from other diagrams [7] . The decay time evolutions for initial 
5° and 5° are [S] 

r I Bl^ J/^/o I = Me-^^' I e^r^*/2(l + COS0,) + e-''^^''^! - cos, 



± sin 0s sin ( Am^ t) J , (1) 

where AF^ is the decay width difference between light and heavy mass eigenstates, AF,, = 
Fl — Fh. The decay width F^ is the average of the widths Fl and Fh, and A/" is a time- 
independent normalization factor. The plus sign in front of the sin 0^ term applies to an 
initial B^ and the minus sign for an initial i?^ meson. The time evolution of the untagged 
rate is then 

F (5° ^ J/V^/o) +F (5° ^ J/V^/o) = Me~^^' I e^^=*/2(l + cos0.) + e-^^=*/2(l-cos 
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Note that there is information in the shape of the hfetime distribution that correlates AF^ 
and (j)s. In this analysis we will use both samples of flavour tagged and untagged decays. 
Both Eqs. [T]and|2]are insensitive to the change 0^ — t- tt — 0s when AF^ — > — AF^. 

2 Selection requirements 

We use a data sample of 0.41 fb~^ collected in 2010 and the first half of 2011 at a centre- 
of-mass energy of 7 TeV. This analysis is restricted to events accepted hy a. J/ip — >■ 
trigger. The LHCb detector and the track reconstruction are described in Ref. [9]. 
The detector elements most important for this analysis are the VELO, a silicon strip 
device that surrounds the pp interaction region, and other tracking devices. Two Ring 
Imaging Cherenkov (RICH) detectors are used to identify charged hadrons, while muons 
are identified using their penetration through iron. 

To be considered a. J/ip ^ f^^f^~ candidate particles of opposite charge are required to 
have transverse momentum, p^, greater than 500 MeV, be identified as muons, and form 
a vertex with fit per number of degrees of freedom (ndof) less than 11. We work in 
units where c = h = 1. Only candidates with dimuon invariant mass between —48 MeV 
to +43 MeV of the J/ip mass peak are selected. Pion candidates are selected if they are 
inconsistent with having been produced at the primary vertex. The impact parameter 
(IP) is the minimum distance of approach of the track with respect to the primary vertex. 
We require that the formed by using the hypothesis that the IP is zero be > 9 for each 
track. For further consideration particles forming di-pion candidates must be positively 
identified in the RICH system, and must have their scalar sum p^ > 900 MeV. 

To select candidates we further require that the two pions form a vertex with a 

< 10, that they form a candidate 5° vertex with the J/ip where the vertex fit x^/ndof 
< 5, that this vertex is > 1.5 mm from the primary, and points to the primary vertex at 
an angle not different from its momentum direction by more than 11.8 mrad. 

The invariant mass of selected n~^fi~7T7r combinations, where the di-muon pair is con- 
strained to have the J/ ip mass, is shown in Fig. |2] for both opposite-sign and like-sign 
di-pion combinations, requiring di-pion invariant masses within 90 MeV of 980 MeV. Here 
like-sign combinations are defined as the sum of 7r+7r+ and tt^tt^ candidates. The signal 
shape, the same for both and _B , is a double-Gaussian, where the core Gaussian's mean 
and width are allowed to vary, and the fraction and width ratio for the second Gaussian 
are fixed to the values obtained in a separate fit to -B^ — ?■ J/ ipcj). The mean values of both 
Gaussians are required to be the same. The combinatoric background is described by an 
exponential function. Other background components are B^ — t- J/iph^ , where h~ can be 
either a K~ or a vr" and an additional 7r+ is found, 5° — )■ J/ipr]', rj' — )■ p7, B^ — )■ J/ip(p, 
— )■ vr+vr"??", and B — )• J/ijjK . The shapes for these background sources are taken 
from Monte Carlo simulation based on PYTHIA fTO] and GEANT-4 [llj with their nor- 
malizations allowed to vary. We performed a simultaneous fit to the opposite-sign and 
like-sign di-pion event distributions. There are 1428±47 signal events within ±20 MeV of 
the B^ mass peak. The background under the peak in this interval is 467±11 events, giv- 
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ing a signal purity of 75%. Importantly, the like-sign di-pion yield at masses higher than 
the gives an excellent description of the shape and level of the background. Simulation 
studies have demonstrated that it also describes the background under the peak. 




5200 5300 5400 5500 



m(J/\j/jt+jt)(MeV) 




5200 5300 5400 5500 

m(J/\j/rt-jt^(MeV) 

Figure 2: (a) Invariant mass of J/ipTi^'K"' combinations when the tt+tt^ pair is required 
to be within ±90MeV of the nominal /o(980) mass. The data have been fitted with 
a double- Gaussian signal and several background functions. The thin (red) solid line 
shows the signal, the long-dashed (brown) line the combinatoric background, the dashed 
(green) line the B~ background (mostly at masses above the signal peak), the dotted 
(blue) line the B — J /ipK background, the dash-dot line (purple) the B — ?► J/ip-K^-K^ 
background, the dotted line (black) the sum of -B^ — )■ J/iprj' and J/ipcj) backgrounds 
(barely visible), and the thick-solid (black) line the total, (b) The mass distribution for 
like-sign candidates. 

The invariant mass of di-pion combinations is shown in Fig. [3] for both opposite-sign 
and like-sign di-pion combinations within ±20 MeV of the candidate mass peak. A 
large signal is present near the nominal /o(980) mass. Other 5° — )• J /%1)tt^tt~ signal events 
are present at higher masses. In what follows we only use events in the /o signal region 
from 890 to 1070 MeV. 
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Figure 3: Invariant mass of vr^vr^ combinations (points) and a fit to the vr^vr^ data 
(dashed hne) for events in the signal region. The region between the vertical arrows 
contains the events selected for further analysis. 

3 S-wave content 



Since the initial isospin of the ss system that produces the two pions is zero, and since 
the G-parity of the two pions is even, only even spin is allowed for the tt'^tt' pair. Since 
no spin-4 resonances have been observed below 2 GeV, the angular distributions are 
described by the coherent combination of spin-0 and spin-2 resonant decays. We use the 
helicity basis and define the decay angles as Oj/^, the angle of the /i"*" in the J/ip rest 
frame with respect to the direction, and 6fg, the angle of the 7r+ in the vr+vr" rest 
frame with respect to the B^ direction. The spin-0 amplitude is labeled as Aqq, the three 
spin-2 amplitudes as A2i, i = —1, 0, 1, and 6 is the strong phase between the A20 ani 
amplitudes. 

After integrating over the angle between the two decay planes the joint angular dis- 
tribution is given by [12] 



00 



dT 



dcos Of^dcosO 



^o + ^^oe''y5(3cos2%„ -1) 



sin^^ 



+ -(|A2i|' + 1^2-1!') (iSsin^^/^cos^^^/J (l + cos^^j/^) 



(3) 



Since the 5° is spinless, when it decays into a spin-1 J /ip and a spin-0 /o, Ojj^ should be 
distributed as sin^6'j/^ and cosOf^ should be uniformly distributed. 

The helicity distributions of the opposite-sign data selected with reconstructed 
J/ipTT^TT' mass within ±20 MeV of the known 5° mass and within ±90 MeV of the nom- 
inal /o(980) mass, are shown in Fig. Ill the data have been background subtracted, using 
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the like-sign data, and acceptance corrected using Monte Carlo simulation. We perform 
a two-dimensional unbinned angular fit. The ratio of rates is found to be 



I A 



20 1 



i^oor 

\A21? + \A2-1? 



(0.1^^:?)%, 
(o.olJ:o')%, 



(4) 



00 1 



where the uncertainties are statistical only. The spin-2 amplitudes are consistent with 
Note that the A20 amplitude corresponds to CP odd final states, and thus would 



zero. 



exhibit the same CP violating phase as the J/ipfo final state, while the A2±i amplitude 
can be either CP odd or even. Thus this sample is taken as pure CP odd. 
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Figure 4: Efficiency corrected, background subtracted angular distributions in the vr+vr" 
mass region within ±90 MeV of 980 MeV and within ±20 MeV of the mass for (a) 
cos^Oj/^, and (b) cosOf^. The solid lines show the expectations for a spin-0 object. 



4 Time resolution and acceptance 

The decay time is defined here as t = md ■ 'p/\'p\^ , where m is the reconstructed 
invariant mass, p the momentum and d the flight vector of the candidate 5° from the 
primary to the secondary vertices. If more than one primary vertex is found, the one that 
corresponds to the smallest IP of the 5° candidate is chosen. 

The decay time resolution probability distribution function (PDF) is determined from 
data using J / if) detected without any requirement on detachment from the primary vertex 
(prompt) plus two oppositely charged particles from the primary vertex with the same 
selection criteria as for J/ipfo events, except for the IP requirement. Monte Carlo 
simulation shows that the time resolution PDF is well modelled by these events. Fig. [5] 
shows the t distribution for our J/tpir^TT^ prompt 2011 data sample. To describe the 
background time distribution three components are needed, (i) prompt, (ii) a small long 
lived background (/lli = 2.64 ± 0.10)% modeled by an exponential decay function, and 
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(iii) an even smaller component (/ll2 = 0.46 ± 0.02)% from 6-hadron decay described by 
an additional exponential. Each of these are convolved individually with a triple-Gaussian 
resolution function with common means, whose components are listed in Table [Tj The 
overall equivalent time resolution is at= 38.4 fs. 

The functional form for the time dependence is given by 



N{t) 



,1 — /lli — /ll2) ■ SG + /lli 



exp(-t/ri) (g) 3G 



+/ll2 ■ [l/r2 ■ exp(-t/r2) ® 3G] . 



(5) 



The fractions /lli and /ll2 , and their respective lifetimes ti and T2, are varied in the 
fit. The parameters of the triple- Gaussian time resolution, 3G, are listed in Table [Tj The 
symbol ® indicates a convolution. 




2.5 3 3.5 4 

Decay time (ps) 



Figure 5: Decay time distribution for prompt J/ipir^TT events. The dashed line (red) 
shows the long lived components, while the solid line (blue) shows the total. 



A decay time acceptance is introduced by the triggering and event selection require 
ments. Monte Carlo simulations show that the shape of the decay time acceptance func 
tion is well modelled by 

[a {t - U)f 



A{t) = C- 



(6) 



l + [a(t-to)]" ' 

where C is a normalization constant. Furthermore, the parameter values are found to be 
the same for simulated B — )■ J/il)K events with K — )■ K^tt^, as for — )■ J/ipfQ. 

Fig. |6|a) shows the J /ipK*'^ mass distribution in data with an additional requirement 
that the kaon candidate be positively identified in the RICH system, and that the K~tt^ 
invariant mass be within ±100 MeV of 892 MeV. There are 36881±208 signal events. 
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The sideband subtracted decay time distribution is shown in Fig. |6](b) and fit using the 
above defined acceptance function gives values of a = (1.89 ± 0.07) ps~^, n = 1.84 ± 0.12, 
to = (0.127 ± 0.015) ps , and also a value of the B lifetime of 1.510±0.016 ps, where the 
error is statistical only. This is in good agreement with the PDG average of 1.519±0.007ps 

m- 




Figure 6: Distributions for B — t- J/tpK* events (a) B candidate mass distribution and 
(b) decay time distribution, where the small background has been subtracted using the 
B candidate mass sidebands. 

Another check is provided by a recent CDF lifetime measurement of B^ — )• J/ipfo of 
l-70^g;}^ ± 0.03 ps obtained by fitting the data to a single exponential [B]. Such a fit to 
our data yields 1.68 ± 0.05 ps, where the uncertainty is only statistical. 

5 Fit strategy 

5.1 Likelihood function characterization 

The selected events are used to maximize a likelihood function 

N 

C = l[P{m„t^,q,), (7) 

i 

where rrii is the reconstructed candidate 5° mass, the decay time, and the total 
number of events. The fiavour tag, g,, takes values of +1, —1 and 0, respectively, if the 
signal meson is tagged as B^, or untagged. The likelihood contains three components: 
signal, long-lived (LL) background and short-lived (SL) background. 
For tagged events we have 

+Ni.i.etlP:!^^{m,)P}'\U) + Nsi.etlP^^^{m{)P!\u), (8) 
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where: (i) P!^{mi) and P^^{mi) are the PDFs describing the dependence on recon- 
structed mass rrii for signal and background events; (ii) P^^^it-, qi) is the PDF used to 
describe the signal decay rates for the decay time t^] (iii) Pl^^it^) is the PDF describing 
the long-lived background decay rates, and Pf^{t^) describes the short-lived background, 
both of which do not depend on the tagging; (iv) e**^^ refers to the respective tagging 
efficiencies for signal, long-lived and short-lived backgrounds. 
For untagged events we have 

PK, t,, 0) = N^,{i - ^)p^^^prit,, 0) 

+NLdl - etl)P:^'{m.)Pl^\u) + iVsL(l - etl)P^J:^{m,)P!\u). (9) 

The total yields of the signal and background components are fixed to the number of 
events determined from the fit to the mass distributions (see Sec. [2]). For both, the PDF 
is a product which models the invariant mass distribution and the time-dependent decay 
rates. The 5° mass spectrum is described by a double-Gaussian for the signal and an 
exponential function for the background (see Fig. [2]). From Eqs. [l]and[2| the decay time 
function for the signal is 

i?(t, qi) oc e~'"°* I cosh — + cos 0^ sinh qiD sin 0^ sin(Amst) 1 . (10) 



The probability of a wrong tag, is included in the dilution factor D = (1 — 2uj) (see 



Section 5.2). 

The signal PDF is taken as a product of the decay time function, R{t,qi), convolved 
with the triple Gaussian time resolution function multiplied with the time acceptance 
function found from J/ipK*^ discussed in Section |4| The background decay time PDFs 
are determined using the like-sign vr^vr^ combinations. The time distribution of the like- 
sign background agrees in both yield and shape with the opposite-sign events in the upper 

mass candidate sideband 50— 200MeV above the mass peak. 

The background functions and parameters are listed in Table [Tj The short-lived back- 
ground component results from combining prompt J/ip events with a opposite-sign pion 
pair that is not rejected by our selection requirements. The long-lived part constitutes 
~85% of the background. 
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Table 1: The PDFs for the invariant mass and proper time describing the signal anc 
background. P^^^ refers to the decay time distribution in Eq. [9 and A is given in Eq. [6 
Where two numbers are listed, the first refers to the 2011 data and the second to the 201Ci 
data. If only one number is listed they are the same for both years. The symbol t refers 
to the true time. 



Pn 



Pt 



Si 

Double-Gaussian {2G) 
2G(m;mo,cri,cr2,/2) 
mo= 5366.5(3) MeV 
ai=8.6(3) MeV 
(72=26.8(9) MeV 
/2= 0.14(2) 



gnal 

Pnt, q) = R{i, q) ® 3G{t - i; fi, a{, a* , a*, 
■A{t] a, n, to) 

/X = -0.0021(1) ps, -0.0011(1) ps 
a{ = 0.0300(4) ps, 0.0295(5) ps 
al/al = 1.92(4), 1.88(3) 
al/al = 14.6(10), 14.0(9) 
/* = 0.23(2), 0.27(3) 

= 0.0136(6), 0.0121(7) 
a = 1.89(7) ps-\ n = 1.84(12), to = 0.127(15) ps 



Long-lived background 



Exponential 




® 2G(t - 


- i; At, (^\, o-2> U)] ■ A{t; a, n, to) 






= 






<y\ 


= 0.088 ps 








= 5.94 ps 








= 0.0137 






r^kg = 0.96 ps 






a = 


= 4.44 ps"\ n -- 


= 4.56, to = Ops 



Short-lived background 



Exponential 



2G(t; fi, a\, a\, f^) ■ A{t; a, n, to) 

All parameters are the same as for LL background 
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5.2 Flavour tagging 

Flavour tagging uses decays of the other b hadron in the event, exploiting information from 
several sources including high transverse momentum muons, electrons and kaons, and the 
charge of inclusively reconstructed secondary vertices. The decisions of the four tagging 
algorithms are individually calibrated using — )■ J/ipK^ decays and combined [H]- 
The effective tagging performance is characterized by e^fgD^, where e^fg is the efficiency 
and D the dilution. We use a per-candidate analysis that uses both the information of 
the tag decision and of the predicted mistag probability to classify and assign a weight 
to each event. The PDFs of the predicted mistag are taken from the side-bands for the 
background and side-band subtracted data for the signal. 

The calibration procedure uses a linear dependence between the estimated per event 
mistag probability i] and the actual mistag probability u given by a; = po +Pi " (^7 ~ 
where po and pi are calibration parameters and (77) is the average estimated mistag prob- 
ability as determined from the calibration sample. In the 2011 data po = 0.384 ± 0.003 ± 
0.009, pi = 1.037 ± 0.040 ± 0.070, and {r]) = 0.379, with similar values in the 2010 sam- 
ple. In this paper whenever two errors are given, the first is statistical and the second 
systematic. Systematic uncertainties are evaluated by using different channels to perform 
the calibration including B — > D*~^fi~u, B^ — )■ J/ipK^ separately from B~ — )■ J/ipK' , 
and viewing the dependence on different data taking periods. For our 2011 sample el^^ is 
(25.6±1.3)% providing us with 365±22 tagged signal events. For signal the mean mistag 
fraction, {■)]), is 0.375±0.005, while for background the mean is 0.388±0.006. After sub- 
tracting background using like-sign events, we determine D = 0.289 leading to an eD^ of 
2.1% P]. 

6 Results 

Several parameters are input as Gaussian constraints in the fit. These include the LHCb 
measured value of Arris = (17.63±0.11±0.02) ps~^ [15j, the tagging parameters po and pi, 
and both the decay width given by the J/ip(j) analysis of = (0.657 ± 0.009 ± 0.008) ps~^ 
and AF, = (0.123 ± 0.029 ± 0.011) ps"^ [16]; we also include the correlation of -0.30 
between F^ and AF^j^ The fit has been validated both with samples generated from 
PDFs and with full Monte Carlo simulations. 

Fig. [7] shows the difference of log-likelihood value compared to that at the point with 
the best fit, as a function of (ps- At each (ps value, the likelihood function is maximized 
with respect to all other parameters. The best fit value is (ps = —0.44 ± 0.44 rad. The 
projected decay time distribution is shown in Fig. [8j 

^The final fitted values of these parameters are shifted by less than 2% from their input values. 
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Figure 7: Log-likelihood profile of 0^ for -8° ~^ J/i'fo events. 



1 ' ' ' ' r 




Decay time (ps) 

Figure 8: Decay time distribution from the fit for J/ipfo candidates. The solid line 
shows the results of the fit, the dashed line shows the signal, and the shaded region the 
background. 

7 Systematic uncertainties 

The systematic errors are small compared to the statistical errors. No additional uncer- 
tainty is needed for errors on Am^, F^, AF^ or flavour tagging, since Gaussian constraints 
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are applied in the fit. Other uncertainties associated parameters fixed in the fit are 
evaluated by changing them by ±1 standard deviation from their nominal values and 
determining the change in fit value of 0^- These are listed in Table |2} An additional 
uncertainty is included due to the possible CP even D-wave. This has been measured at 
(0.0]'^o!o)% of S-wave and contributes a small error to 0s, +0.007rad, as determined by 
repeating the fit with the mistag rate increased by 1.7%. The asymmetry in production 
between 5° and 5° is believed to be small, about 1%, and similar to the same asymmetry 
in production which has been measured by LHCb to be about 1% [IT]. The effect 
of neglecting a 1% production asymmetry is the same as ignoring a 1% difference in the 
mistag rate and causes negligible bias in 0^. 

Table 2: Summary of systematic uncertainties. Here A^bkg refers to the number of back- 
ground events, iVsig the number of signal, A^,;' the number of 77', a the exponential back- 
ground parameter for the Bl candidate mass, iVLL/Abkg the long-lived background frac- 
tion. The Gaussian signal parameters are the mean mo, the width cr(m); and n 
are the three parameters in the acceptance time function. The resolution in signal time 
is given by cr(t), and the background hfetime by rbkg. The final uncertainty is found by 
adding all the sources in quadrature. 



Quantity (Q) 


±AQ 


-|- Change 


—Change 






in 0s 


in (ps 


A'bkg 


10.1 


0.0025 


-0.0030 




3.4 


-0.0001 


-0.0001 


A^sig 


46.47 


-0.0030 


0.0028 


a 


1.7-10"^ 


-0.0002 


-0.0002 




0.0238 


0.0060 


-0.0063 


mo (MeV) 


0.32 


-0.0003 


0.0011 


a{m) (MeV) 


0.31 


-0.0026 


0.0020 


Thkg (ps) 


0.05 


-0.0075 


0.0087 


a{t) (ps) 


5% 


-0.0024 


0.0022 


to (ps) 


0.015 


0.0060 


0.0050 


a (ps~^) 


0.07 


-0.0065 


-0.0065 


n 


0.12 


-0.0089 


-0.0089 


CP-even D-wave 




0.0070 





Total Systematic Error 


+0.018 


-0.017 



8 Conclusions 

Using 0.41 fb^^ of data collected with the LHCb detector, the decay mode — )■ J/ipfo, 
fo — 7- 7r+7r~ is selected and then used to measure the CP violating phase, 0^. We perform 
a time dependent fit of the data with the 5° lifetime and the difference in widths of the 
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heavy and light eigenstates constrained. Based on the hkehhood curve in Fig. [7] we find 

(p, = -0.44 ± 0.44 ± 0.02 rad, 

consistent with the SM value of — 0.03631q oqi5 rad [IJ. Assuming the SM , the probability 
to observe our measured value is 36%. There is an ambiguous solution with (j)s ^ tt — 4>s 
and ATs — )■ — AF^. The precision of the result mostly results from using the tagged 
sample, though the untagged events also contribute. 

LHCb provides an independent measurement of 0s = 0.15 ± 0.18 ± 0.06 [T6] using the 
— )■ J/ip(j) decay. Combining these two results, taking into account all correlations by 
performing a joint fit, we obtain 

(j)s = 0.07 ± 0.17 ± 0.06 rad (combined). 

This is the most accurate determination of 0^ to date, and is consistent with the SM 
prediction. 
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